Effect of controlled corrugation on capillary condensation of colloid-polymer mixtures 



Andrea FortinQ and Matthias Schmidt 
Theoretische Physik II, Physikahsches Institut, Universitat Bayreuth, 
Umversttatsstrafle 30, D-954i^ Bayreuth, Germany 

We investigate with Monte Carlo computer simulations the capillary phase behaviour of model 
colloid-polymer mixtures confined between a flat wall and a corrugated wall. The corrugation is 
modelled via a sine wave as a function of one of the in-plane coordinates leading to a depletion 
attraction between colloids and the corrugated wall that is curvature dependent. We find that for 
increased amplitude of corrugation the region of the phase diagram where capillary condensation 
occurs becomes larger. We derive a Kelvin equation for this system and compare its predictions 
to the simulation results. We find good agreement between theory and simulation indicating that 
the primary reason for the stronger capillary condensation is an increased contact area between 
the fiuid and the corrugated substrate. On the other hand, the colloid adsorption curves at colloid 
gas-liquid coexistence show that the increased area is not solely responsible for the stronger capillary 
condensation. Additionally, we analyse the dimensional crossover from a quasi-2D to a quasi- ID 
system and find that the transition is characterised by the appearance of a metastable phase. 



I. INTRODUCTION 

The equilibrium behaviour of a fluid in contact with 
a solid substrate is ruled by the interfacial free energy, 
i.e., the amount of free energy needed to create the in- 
terface [1]. Since this free energy is the product of the 
interfacial tension and the total contact area, the fluid 
equilibrium behaviour at a substrate can be controlled 
by either manipulation of the chemical properties of the 
surface (interfacial tension) or the geometry of the sur- 
face (contact area) . It is known that corrugated surfaces 
have a higher contact area with respect to a flat surface 
with the same cross-sectional area and in nature, 

corrugation gives the Lotus flower ^ its characteristic 
hydrophobicity. Efforts to mimic the Lotus effect resulted 
in the production of micromachined surfaces whose wet- 
tability has been controlled by proper surface microstruc- 
turing O Ej . From a fundamental point of view wetting 
and capillary condensation on structured [7] or curved 
substrates [SI 12] as well as on surfaces with wedge geome- 
try dnmi] have been studied in detail (see also the review 
by Bonn et al. [12] and references therein). Nevertheless, 
while a lot of research has been done for simple fluids, 
the effect of corrugation has been neglected in complex 
colloidal fluids because the roughness of a substrate nor- 
mally occurs at length scales that are much smaller than 
the typical size of a colloidal particle. 

Recently, a technique has been introduced that allows 
the controlled wrinkling of surfaces [Ol [H] on the mi- 
cron scale. The effect of wrinkled substrates on crys- 
tallisation has been investigated [THj, and the resulting 
structures have been used to enhance Surface Raman 
Spectroscopy [TS] . Clearly, the wrinkling technique could 
allow a systematic study of the effect of roughness on the 
phase and wetting behaviour of complex fluids. 

In this article, we investigate with computer simula- 
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tions the phase behaviour of a mixture of colloid and non- 
adsorbing polymers conflned between a wrinkled wall and 
a flat wall. We chose colloid-polymer mixtures because 
they have been extensively studied in the literature [171 - 
both from a theoretical and experimental point of 
view. In particular, such mixtures allowed the investiga- 
tion of a variety of fundamental inhomogeneous statisti- 
cal phenomena such as wetting at single flat walls [5D1 - 
I34| and capillary condensation and evaporation between 
flat walls [35ll40j . Here, we describe the colloid-polymer 
mixture using the Asakura-Osaawa-Vrij (AOV) [U 22] 
model. Interestingly, the depletion attraction between 
colloids and the corrugated wall depends on the local 
curvature of the wall (the same effect occurs for mix- 
tures of hard spheres [331 SI])- We flnd that for in- 
creased amplitude of corrugation the capillary conden- 
sation is stronger, i.e. the region of the phase diagram 
(phase space) where capillary condensation occurs be- 
comes larger. We also analyse the conflned system us- 
ing a simple thermodynamic description, derive a Kelvin 
equation and study its validity by comparing the theoret- 
ical predictions to the simulation results. The agreement 
between theory and simulation is good, and indicate that 
the increased contact area between the fluid and the sub- 
strate (due to corrugation) is the primary cause for the 
stronger capillary condensation. A dramatic increase in 
colloidal adsorption at the corrugated wall is also primar- 
ily due to the increased surface area. On the other hand, 
further analysis suggests that other effects influence the 
thermodynamics of the system, like for example the cur- 
vature dependence of the particle wall interaction. Fur- 
thermore, we flnd that the crossover from a quasi-2D to 
a quasi- ID system leads to the emergence of metastable 
phases characterised by a sequence of fllled and empty 
quasi-lD channels. 

The paper is organised as follow. In Sec. |Tl] we intro- 
duce the model and the simulation technique. In Sec. IV 
we derive an expression for the Kelvin equation for cor- 
rugated walls, while in Sec. |V] we present the simulation 
results. Finally in Sec. |VI| we draw our conclusions. 



2 




Figure 1. a) Schematic drawing of the simulation box and model. Colloids of diameter (Tc and spheres with diameter ap 
representing the polymers are confined between hard walls: a flat wall at the bottom [z — 0) with area Aflat = x Ly, and 
a sinusoidal wall of wavelength A and amplitude a at the top {z = H + asin(27rx/A)) . b) Cross section of the simulation box. 
The depletion zones at the walls are delimited by dashed lines. The grey area indicates the volume accessible to the particles, 
c) Colloids and Polymer at the flat wall. The depletion zones are delimited by dashed lines. The grey area indicates the gain 
in free volume due to the overlap of two depletion regions, d) Like c) but for a colloid at a sinusoidal wall. The gain in free 
volume depends on the local curvature of the wall. 



II. MODEL 

The Asakura-Osaawa-Vrij (AOV) [HI |32] model is a 
binary mixture of colloidal hard spheres of diameter CTc 
and of spheres with diameter dp representing polymer 
coils (Fig. [T^). The pair interaction between colloids is 
that of hard spheres: Vcdr) = cx) if r < ctj. and zero 
otherwise, where r is the separation distance between 
particle centres. The interaction between a colloid and 
a polymer is also that of hard spheres: Vcp{r) = cxi if 
r < (ac + <Jp)/2 and zero otherwise. The polymers, how- 
ever, are assumed to be ideal, hence the polymer-polymer 
interaction vanishes for all distances, Vpp{r) = 0. The 
size ratio q = ap/cTc = 1 is a geometric control parame- 
ter. 

The colloid-polymer mixture is confined between one 
smooth planar hard wall at z = and one sinusoidal hard 
wall ed, z = H + asin(27ra;/A), where H is a measure of 
the distance of the two walls, a is the amplitude of the 
sinusoidal corrugation, A is its wavelength, and x is the 
lateral coordinate perpendicular to the wrinkles. Any 
overlap between particles and walls is omitted, therefore 
the external potentials acting on species i = c,p 



ri e Vhcc 
oo otherwise, 



(1) 



where Vfroo is the free volume available to the particles 
as sketched in Fig.[I]D). 

The polymers induce an effective colloid-colloid and 
colloid- wall attraction that is of entropic origin and is due 



to the so-called depletion effect. In Fig. [T]:), we show a 
colloid in contact with a wall. Around the colloid there is 
a depletion region (dashed circle) prohibited to the poly- 
mers due to the colloid-polymer hard-core interaction 
(the chains cannot penetrate the colloids) . The hard-core 
interaction between polymers and the wall gives rise to a 
depletion zone at the wall (dashed line in Fig. [ij;) . If two 
colloids approach each other, so that two depletion zones 
overlap there is an increase in free volume for the poly- 
mer chains, i.e., an increase in entropy [T71 [THl HH HI]- 
The increase in entropy can be described by an attractive 
interaction between colloidal particles. Likewise, there is 
an increase in entropy when the depletion zones of col- 
loids and wall overlap (grey zone). The larger the ex- 
cluded volumes, the stronger is the effective attraction 
between the hard wall and the colloids. Figure [iji) il- 
lustrates the depletion zone of a curved wall. Clearly, 
the gained free volume (grey zones) depends on the local 
curvature of the wall. The concave part of the wall leads 
to a larger gain in free volume (more attractive) than the 
convex part of the wall (less attractive). 

We denote the packing fractions by rji = 
■KcrfNi/ [QAH)^ where Ni is the number of particles 
of species i and A is the cross-sectional (normal to the 
z-direction) area of the wrinkled wall, and also equal 
to the area of the flat wall. Additional thermodynamic 
quantities are the scaled chemical potential /3/ic and /3/ip 
of colloids and polymers, respectively, and the polymer 
reservoir packing fraction ff!^ = exp(/3/ip) of a reservoir of 
pure polymers that is in chemical equilibrium with the 
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system. The inverse temperature is /3 = l/ksT, with ks 
the Boltzmann constant and T the temperature 



III. SIMULATION METHOD 

Fig. [T^) displays an illustration of the model, for which 
we carried out Monte Carlo computer simulations for a 
wide range of different values of particle concentrations 
and for several values of the amplitude a. The wavelength 
of the corrugation was fixed to the value of A = lOtXc. Pe- 
riodic boundary conditions were applied in the x- and y- 
directions, and the box size in the x-direction was chosen 
as 4A. 

We carried out Monte Carlo simulations in the grand 
canonical ensemble, i.e. with fixed volume, temperature, 
and chemical potentials /3^c and /3^p . 

To study the phase coexistence, we sample the proba- 
bility P{Nc)\f_i^^f_i^ of observing Nc colloids in a volume V 
using the successive umbrella sampling [45]. We use the 
histogram reweighing technique to obtain the probability 
distribution for any fi'^ once P{Nc)\^i^.^j, is known for a 
given 

lnP{N,)\^,^^^ = \nP{N,)l^,^^ + NM - Pt^c). (2) 

At phase coexistence, the distribution function P{Nc) 
becomes bimodal with two separate peaks of equal area 
for the colloidal liquid and gas phases. We determine 
which fi'^ satisfies the equal area rule 



P{N,)\^,,,,^dN,^ P{N,)l,,.,,^dN,, (3) 

-'(Na) 



with the average number of colloids 



N,P{N,)\^.dN, 



(4) 



using the histogram reweighing equation ([2]). The sam- 
pling of the probability ratio P{N)/ P{N + 1) is done, in 
each window, until the difference between two successive 
samplings of the probability ratio is smaller than Ix 10"''. 
To improve the sampling accuracy, we used the cluster 
move introduced by Vink and Horbach [46) . 



IV. THEORY: KELVIN EQUATION 

We apply a thermodynamic treatment in the limits of 
A/(T ^ 1 and H/a ^ 1 to the system sketched in Fig. n 
The derivation follows closely the derivation of Evans [47 
for fluids between smooth parallel walls. For simplicity 
we limit the theoretical derivation to the case of one- 
component fluids. As will become clear at the end the 
section, the result can be generalised effortlessly to multi- 
components fluids. 



The free energy in the grand canonical ensemble is the 
grand potential fl. For a fluid between a top and a bot- 
tom walls of area Ai and A2, respectively, the grand po- 
tential is the sum of the bulk free energy per unit vol- 
ume Lu and the surface energy contributions from the two 
walls. 



(5) 



where V is the volume, fi is the chemical potential and 
7'^,7^, are the fluid- wall interfacial tensions for the top 
and bottom wall, respectively. 

We introduce the bulk chemical potential at coexis- 
tence fj,b, and use it as a reference state for the chemical 
potential of the confined system 

Assuming A/i small we Taylor expand equation ([S]) 
around /Zb, recalling that the bulk density is 

dui 

^ djjL 



and that the adsorption is 
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The Taylor expansion yields 



V 



Ai 



V 



V 



2„2 



V 



7' (Mb) 



r2(^,)A^.(6) 



The coexistence between a liquid and a gas phase in- 
side the slit occurs when both thermal and mechanical 
equilibrium conditions are satisfied, i.e., the chemical po- 
tential and pressure of the gas and liquid phases are the 
same. These two conditions are satisfied when the grand 
potential of the two phases is the same — r2'"i(/i). 

Given that = uj^^^{i.ib), the above equilibrium 

condition inside the slit leads to the following relation 



{pUq - Pgas)A/X + 

V 

^2 /2 2 



^ \ /gas /liq/ ' V gas ^ liq/ 



(r^ _ r2 )Am 



(7) 



We stress that the above relation ^ is valid for any type 
of surface, as we only requested that A/z be small. 

The previous relation greatly simplifies if we further 
assume that the adsorption at the wall is negligible and 
that 7"^ = 7^ = 7, in which case 



Am = ( 



Al ^ ^2 J7liq - 7gas) 
V V (pliq - pgas) 



(8) 



We stress that the assumption 7"^ = 7^ = 7 implic- 
itly considers a curvature independent interfacial tension. 
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Figure 2. The roughness ratio R for a sinusoidal wall, a) For 
fixed amplitude a = 5(Tc and changing wavelength, b) For 
wavelength A = lOtJc and changing amplitude. 



That is, we treat the limit of a/ a ^ 0. This appar- 
ently crude approximation is justified in the framework 
of the Kelvin equation that treats systems in the limit 
H/a ^ 1. Strictly speaking, the Kelvin equation is 
not valid for strongly confined systems, but has been 
shown before that its predictions remain qualitatively 
good down to quasi-2D systems. This theoretical descrip- 
tion would on the other hand break down in the quasi- ID 
limit, because it would predict a phase transition that in 
reality does not exists for a bulk one-dimensional system. 



A. Two flat walls 

If we consider two flat walls of equal area Ai = A2 — 
A^''\ and volume V = A^^^H we obtain the standard 
equation [47] 



2 (7liq-7gas) 
H (piiq - Pgas) 



(9) 



B. One corrugated wall and one flat wall 

A system with one flat wall and one corrugated wall 
has a volume V = A^^*^{H), with {H) the average sepa- 
ration distance between the walls. The relation (|8| then 
becomes 



corr /I 



R) 



(7liq - 7gas) 
(Pliq - Pgas) ' 



(10) 



R = 



is the ratio between the flat and corrugated area and is 
equivalent to the Wenzel roughness ratio defined as the 
ratio between the contact area and the geometric cross 
sectional area [2]. 

By comparing Eqs. ([9| and ( 10 ) we find also that the 
ratio between the chemical potential shifts of the corru- 
gated and flat walls is related to the Wenzel ratio via 



1 + R 



(11) 



C. Two corrugated walls 

For the case of two corrugated walls the general Kelvin 
equation ([8| becomes 



^^corr _ 2 ^ (7liq 7gas) 
{H) (Pliq - Pgas) ' 



(12) 



By comparing Eqs. ( 12 1 and ([£]) we find that for the case 
of two corrugated walls the ratio is 



A^™"7A^"''* = R 



D. Generalisation to binary mixtures 

The derivation of the Kelvin equation for binary mix- 
tures [3S] proceeds along the lines described above for 
a one-component fluid. A closed set of equation for 
the chemical potential shifts A/Xc, and Afj,p, for colloid 
and polymer, respectively, can however be obtained only 
when an independent relationship between A^c, and Afip 
is used. In Ref. l37l three possible choices for the rela- 
tionship are outlined. Independent of the choice of the 
relationship, the equations for binary mixtures for corru- 
gated walls leads to 



A/i™"/A/.; 



flat 



flat 



1 + R 
2 

1 + R 



(13) 



and 



A^:;°'VA^;^^* = R 



(14) 



A^i^-VAAi^^* - R , 
for one and two corrugated walls, respectively. 
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E. Roughness ratio for a sinusoidal wall 

We consider a wall corrugated in one direction by 
a sinusoidal functional form as shown in the sketch of 
Fig. The cross sectional area is j4flat = L^r x Ly, 

where and Ly are the parallel and perpendicular di- 
rections with respect to the sinusoidal direction, and 
L-c = nX, where n is an arbitrary number and A is the 
wavelength. Therefore, the roughness ratio can be writ- 
ten as 



R = 



1 



^COs2(?f.x) 



(15) 



where the numerator is the line integral over the sinu- 
soidal path. The elliptical integral can be solved numer- 
ically. The value of the integral ( 15 1 is the same for any 



integer value of n, but changes when non-integer values 
of n are chosen. Therefore, an explicit dependence of the 
integral on n is left for the general case. 

Figure [2] shows the behaviour of the roughness ratio 
defined in Eq. ( |15[). For fixed amplitude and increasing 
wave length (Fig. |2k) the ratio R decreases. Clearly for 
A — cx) the ratio goes to one. On the other hand, for fixed 
wavelength A (Fig. [2]d) the roughness ratio R increases 
monotonically for increasing values of the amplitude of 
corrugation a. 
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Figure 3. Phase diagram in the {/Sfip, Pfic) representation. 
Shown are the binodals for wavelength A = 10a c and ampli- 
tudes a/ac~0, 5, 7, 9, 11, 13, and 15. Also shown are the 
binodals for the bulk system (dashed line) and for the system 
confined between two fiat walls at distance H = 5(Tc (thick 
continuos line). 



V. SIMULATION RESULTS 
A. Phase diagram 

Figure [S] shows the phase diagram in the (/3^p, (3fJ-c) 
representation. Each line represent a binodal, i.e. the 
value of colloid and polymer chemical potentials at which 
phase separations occurs. For chemical potentials in the 
region above a binodal the gas phase is stable, whereas 
in the region below the binodal the liquid phase is sta- 
ble. In particular, the capillary binodals for corrugation 
amplitude a/a,, = 0—15 are shown (symbols). As refer- 
ence we also show the bulk binodal (dashed line) , and the 
capillary binodal for flat walls with separation distance 
H/(Tc = 5 (thick continuous line). The binodal for the 
colloid-polymer mixtures confined between two flat walls 
(i.e., a = H/ac = 15) is shifted toward smaller col- 
loidal chemical potentials and higher polymer chemical 
potential with respect to the bulk binodal (thin dashed 
line). The shift of the binodals indicate the occurrence 
of capillary condensation, because in the region between 
the confined and bulk binodals, the stable phase for the 
confined system is the liquid, while for the bulk system is 
the gas. This effect was extensively studied in the litera- 
ture [351 133 IMl EH] for colloid-polymer mixtures confined 
between flat walls. 

Here, for colloid-polymer mixtures confined between 
one flat and one corrugated wall, we find that for in- 
creasing corrugation, i.e., increasing amplitude a of the 



sinusoidal wall, the binodals shifts toward larger polymer 
chemical potentials and smaller colloid chemical poten- 
tials, indicating progressively stronger capillary conden- 
sation with respect to the flat walls, i.e., larger region 
of phase space where the capillary condensation occurs. 
The critical points were not calculated because an ac- 
curate determination would require a careful flnite-size 
analysis [38] that is beyond the scope of the current work. 
The end-points in Fig. [3] were determined as the lowest 
value of the chemical potentials at which a double peak 
in the probability distribution was observed. 

One trivial interpretation for the behaviour is that the 
sinusoidal wall introduces a length scale in the system 
that is smaller than the average wall separation distance, 
namely /imm = H — a, that strengthens the capillary con- 
densation. However, this interpretation is not supported 
by our results. For example for a =^ Hctc, the smallest 
length scale is hmin — 15(Tc — Hcc — 4(Tc- As a refer- 
ence in Fig. [3] we plot the binodal for two flat walls at 
distance H/ac — 5. The a = 11(7^ binodal (left triangles) 
is clearly separated from it. 

For amplitude a — 15(Jc the system consist of a se- 
ries of independent quasi-lD channels. Interestingly, we 
find that the binodal lines for increasing amplitudes a, 
slowly approach the binodal of the system with a single 
groove, therefore the system undergoes a slow dimen- 
sional crossover from a quasi-2D to a quasi-lD system. 
We stress that in one-dimensional systems the gas-liquid 
phase transition does not exist in the thermodynamic 
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Figure 4. The ratio A^c°''7A^ic * for different amplitudes 
a of the sinusoidal wall at fixed wavelength A = 10a c- We 
compare results of the computer simulations (symbols) for 
different values of polymer reservoir packing fraction rfp with 
the theoretical prediction (line) of equation (10 1. 



limit. Even in quasi-lD pores a rounding of the tran- 
sition is noted jSOj and a multi domain structure is ob- 
served instead of a proper gas-liquid separation. 



B. Comparison with the Kelvin equation 



The Kelvin equation (11) suggests, on the other hand 
that the stronger capillary condensation for increasing 
amplitudes a is due to the increased surface area of the 
corrugated wall. To evaluate the validity of this inter- 
pretation we compare theoretical results with the sim- 
ulation findings. In particular, we compare the ratio 
A/z™''''/A/i^^' found in simulations to the theoretical re- 
sult of Eq. ( |Tl| ). Remarkably, the simulation results 
start at smaller values of the ratio but they follow the 
same trend as the theory results for increasing amplitudes 
a/ac, demonstrating that the stronger capillary conden- 
sation is mainly driven by the increased surface area of 
the corrugated wall. We believe that the smaller values 
of the ratio A/x™''''/A/i^'** in simulations with respect to 
the theory are due to the large adsorption of colloids at 
the sinusoidal wall in the gas phase. The Kelvin equa- 
tion also neglects the curvature. Recent works on curved 
surfaces [53 suggest that the inclusion of curvature 
effects would lead to smaller values of the 7iiq — 7gas dif- 
ference. This leads to a decrease of the shift in chemical 
potential predicted by the Kelvin equation, leading to a 
better agreement with simulation results. 



C. Adsorption 

Further insights can be gained by calculating the col- 
loid adsorption F for state points at gas coexistence. First 
of all, we carried out simulations of colloid-polymer mix- 
tures between two flat walls at separation H — 15ctc to 
measure the adsorption at one flat wall 



Ffla 



N, 



flat 



N, 



bulk 



2^flat 



(16) 



at the gas coexistence state points of the phase diagram 
shown in Fig. [3] Subsequently, we performed simulations 
for the mixtures confined between one flat and one si- 
nusoidal wall and computed the adsorption at a single 
sinusoidal wall 



F. 



N, 



bulk 



^fla 



Ffla 



(17) 



Here we choose as reference the cross-sectional area 

^flat 

as it is common for corrugated surfaces where the real 
contact area is not easily known. We find that the adsorp- 
tion increases dramatically with increasing amplitudes a 
and changes little for increasing polymer reservoir pack- 
ing fractions (Fig. [5^) . 

In order to evaluate the effect of the chosen reference 
area, we also calculated the adsorption using as a refer- 
ence the real contact area A^°^^ , i. e. 



corr 



^corr _ ^bull 



A'' 



Ffla 



(18) 



The increase in colloid adsorption for increasing ampli- 
tudes a is now less dramatic (Fig. ^Ip) , note the different 
y-axes scale) but still clearly visible. Therefore, the in- 
creased surface area of the sinusoidal wall is not the only 
factor responsible for the increased adsorption and con- 
sequently for the stronger capillary condensation. 

Interestingly, the simulation snapshots for state points 
at gas coexistence for polymer reservoir packing fraction 
rjp — lA (shown in Fig. 6]) show that colloids are adsorbed 
inside the grooves of the sinusoidal walls, that is in the 
concave parts of the corrugated wall. 



D. Free energy curves 

We next explore the dimensional crossover that oc- 
curs at amplitudes a/dc comparable to the wall dis- 
tance H. Consider that the free energy (grand poten- 
tial) difference with respect to a reference state is given 
by -ln(P(7/c)Uc,A<p): where P(?7c)U, ,fj.p) i*^ the probabil- 
ity of observing the system at colloid packing fraction ijc- 
In Fig. [T^) we plot the free energy curves for amplitudes 
a=l, 3, 5, 7 (Tc. The two minima correspond to the gas 
and liquid coexisting densities. For intermediate pack- 
ing fractions rjc the free energy decreases for increasing 
packing fraction, as expected for systems in slit geom- 
etry. In Fig. [7]d) the free energy curves for amplitudes 
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Figure 5. Adsorption at the corrugated wall for increasing 
values of the corrugation amplitude a/oc- a) Colloidal ad- 
sorption Fcorr with the cross-sectional area A^'^^ as reference 
surface. (Eq. (|17[)). b) Colloidal adsorption T'o^r with the 
real contact area ^1'^°'^'^ as reference surface (Eq. (18l). 



a=9, 11, 13, 15 (Tc are shown. The approaching dimen- 
sional crossover is characterised by the appearance of a 
sequence of minima in the free energy curves. In our 
system the local minima appear first for o/ctc = 9 and 
becomes stronger and stronger for increasing amplitudes. 
These minima indicate the presence of metastable phases 
characterised by coexistence of channels filled with fluid 
phase and channels containing the gas phase. Each time 
a channel is completely filled we find an energy minimum 
due to the decreased interfacial energy as shown in Fig. [H] 
for polymer reservoir packing fraction rfp — 1.4 and cor- 
rugation amplitude a/uc = 11. These phases are similar 
to the zehra phase predicted for optically confined mix- 





Figure 6. Snapshots of the simulation at gas-liquid coex- 
istence for the packing fraction corresponding to the equi- 
librium gas phase and polymer reservoir packing fraction 
rip = 1.4. The colloidal particles (cyan/light grey spheres) 
and the polymers (green/dark grey spheres) are confined be- 
tween one flat wall (bottom) and one sinusoidal wall (top). 
The walls occupy the white sections of the figures, a) Ampli- 
tude a/(Tc = 1 b) Amplitude a/(Tc = 5 c) Amplitude a/ac = 9 
d) Amplitude a/ac = 13 



tures [53H55 



VI. SUMMARY AND CONCLUSIONS 

We traced the phase diagram of model colloid-polymer 
mixtures confined between one fiat wall and one corru- 
gated wall with grand-canonical MC simulations. The 
corrugation was modelled by a sinusoidal function of am- 
plitude a. We found that for increasing values of a the 
capillary condensation get stronger, i.e., the region of 
parameter space where the capillary condensation oc- 
curs becomes larger. We derive a Kelvin equation for 
the system that predicts that capillary condensation in 
the system is controlled by the Wenzel ratio R — 
between areas of the corrugated and flat walls. We 
very good agreement between simulation results and the 
Kelvin equation prediction, indicating that the increased 
contact area between the fluid and the substrate (due to 
corrugation) is the primary cause for the stronger capil- 
lary condensation. The analysis of the adsorption of col- 
loidal particles shows a strong preference of the colloids 
to adsorb at the corrugated wall in agreement with the 
stronger capillary condensation. The finding is also cor- 
roborated by a visual inspection of the simulation snap- 
shots, which show a strong adsorption of colloids in the 
grooves of the corrugated walls. Nevertheless, we find 
an increase of adsorption for increased surface area also 
when the real contact area is used as reference. One pos- 
sible explanation for this effect is that the depletion effect 
responsible for the colloid-wall effective attractive inter- 
action is curvature dependent. The overlap area between 
the particle and wall excluded volumes is different for 
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teraction. We also observe a dimensional crossover from 
a quasi-2D system to a quasi-lD system. The quasi-lD 
system occurs when the amplitude of the corrugated wall 
is equal to the average separation distance, i.e., all chan- 
nels are decoupled. The free energy curves show that the 
approaching crossover is characterised by the appearance 
of a metastable phase, with partially filled channels, sim- 
ilar to the 'zebra' phase predicted for optically confined 
colloid-polymer mixtures |53l I54j . 

In this paper we only considered capillary condensation 
a) ■ b) ■ 




Figure 8. Snapshots of the simulation at polymer reservoir 
packing fraction rjp = 1.4 and amplitude a/ac = 11. The 
colloidal particles (cyan/light grey spheres) and the polymers 
(green/dark grey spheres) are confined between one flat wafl 
(bottom) and one sinusoidal wall (top). The walls occupy 
the white sections of the figures, a) One channel filled with 
colloids at state point a in Fig. [Tji). b) Two channels filled 
with colloids at state point P in Fig. [tJj). c) Three channels 
filled with colloids at state point 7 in Fig.jTjs). d) Equilibrium 
liquid phase. 



Figure 7. Logarithm of the probability distribution (free en- 
ergy) at polymer reservoir packing fraction rjp = 1.4 and col- 
loid chemical potentials at coexistence, a) For amplitudes 
a=l, 3, 5, 7 (Tc. b) For amplitudes a— 9, 11, 13, 15 CTc. 

different wall curvatures. Therefore the wall is more at- 
tractive at the concave side and less attractive at the con- 
vex edge. In our simulation we used a fixed wavelength 
and changing amplitudes, i.e., the curvature of the corru- 
gated walls would be different for the different amplitudes 
leading to a different effective wall-colloid attractive in- 



and neglected two other relevant phenomena, namely 
wetting [12] and wedge filling transitions [TTl [SHI 157] . 
Given the finite system size as well as the finite wedge 
geometry our assumption is completely reasonable, but 
an extension of this work would be to analyse the wet- 
ting and wedge filling behaviour of a single corrugated 
wall and compare to works on wedge filling transition in 
finite geometries [58H60] . 

A recent paper [55] on fiuids in periodic confinement 
showed unexpected correlations between the channels 
that deserve further investigation in the current system. 
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